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Abstract 

In this work we address the "smoking ring" propulsion technique, originally proposed 
by Purcell in [1] . We first consider self-locomotion of a doughnut-shaped swimmer powered 
by surface tank-treading. Different modes of surface motion are assumed and propulsion 
velocity and swimming efficiency are determined. The swimmer is propelled against the 
direction of its outer surface motion, the inner surface having very little affect. The simplest 
swimming mode corresponding to constant angular velocity, can achieve propulsion speeds 
of up to 66% of the surface tank-treading velocity and swimming efficiency of up to 13%. 
Higher efficiency is possible for more complicated modes powered by twirling of extensible 
surface. A potential practical design of a swimmer motivated by Purcell's idea is proposed 
and demonstrated numerically. Lastly, the explicit solution is found for a two-dimensional 
swimmer composed of two counter-rotating disks, using complex variable techniques. 

1 Introduction 

Nano-technology is about the control of tiny objects. The benefits of tiny artificial swimmers 
for medicine, for example, could be enormous [2j. A tiny robot may swim through the arteries, 
digestive system, spinal canal, etc., transmit images and deliver microscopic payloads to parts of 
the body, or perform some kind of therapeutic action outside the reach of existing technologies. 

Tiny swimmers, be they micro-organisms or microrobots live in a world dominated by fric- 
tion pp. Swimming at low Reynolds numbers is associated with energy dissipation and because 
of the absence of inertia, is often counter-intuitive. Micro-organisms had eons of evolution to 
optimize their locomotion and so are a natural source of inspiration and imitation in design of 
artificial microswimmers (e.g. [3])- Nevertheless, it is natural to ask whether the mechanisms 
used by micro-organisms are also the method of choice for artificial tiny devices. Microor- 
ganisms are different from artificial devices both in structure, time scales and functions, so a 
mechanism that is optimal for one may be inappropriate for the other. For example, useful 
robots would need to swim much faster than microorganisms. This implies that robots would 
be energy guzzlers compared with organisms which enjoy the additional luxury of swimming in 
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Figure 1: Low- Reynolds-number-swimmers proposed by Purcell: (A) "smoking ring" propulsion 
of a torus; (B) two cells stuck together and rolling on one another (redrawn from pQ) Arrows 
pointing to the left show the direction of propulsion, while smaller arrows at the swimmers' 
surface show the direction of the tank-treading motion. 

a sea of nutrients. Autonomous artificial swimmers need also carry on "battery" and/or "en- 
gine" , therefore, the issue of swimming efficiency becomes central. Swimming or hydrodynamic 
efficiency rates different swimmers in terms of the power invested in swimming. 

Swimming of micro-organisms has a rich history in applied mathematics [3] and several 
modes of natural locomotion, such as "flexible oar" (i.e. beating elastic tail) 0E1I7], the "cork- 
screw" (i.e. rotating helical flagellum) [21 El El El E] , cilia [T2] or surface waves [131 E] are 
quite well understood. Artificial biomimetic (e.g. nature-inspired) propellers, powered by either 
beating or rotating tails/filaments [151 El El El E] were recently fabricated and tested. As 
the above natural propulsion techniques are characterized by quite low hydrodynamic efficiency, 
it is both conceptually interesting, and technologically important to understand what strategies 
may lead to effective propulsion in a setting dominated by viscous friction. Since tiny swimmers 
are a challenge to build while the theoretical framework is well understood it makes sense to 
put some effort into toying with mathematical models even if, at the moment, there may be 
no clear concept for their realization. Theoretical works on propulsion of artificial swimmers, 
that are not necessarily inspired by natural modes of locomotion, at low Reynolds numbers is 
a young enterprize and some mechanisms such as three-link Purcell's swimmer [201 [21] three- 
linked-sphere swimmer [22], swimmer propelled by arbitrary non-retractable periodic shape 
strokes [Ml 123], "pushmepullyou" [25] and others were recently studied. 

A particularly promising class of strategies is where the motion is, in a sense, only apparent; 
where a shape moves with little or no motion of material particles. In our recent work [26] we 
focused on locomotion powered by surface treadmilling. In surface treadmilling the elongated 
microswimmer with length £ and thickness d moves without a change of shape, by a tangential 
surface motion. Surface is generated at the front end and is consumed at the rear end. In 
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contrast to actin and microtubules, the surface treadmilling does not rely on the exchange 
of material with the ambient fluid. (The swimmer needs, of course, an inner mechanism to 
transfer material from its rear to its front which we did not account for in the model.) As the 
slenderness increases, the hydrodynamic disturbance created by the surface motion diminishes, 
i.e. the swimmer is propelled forward with the speed of boundary displacement backward, as the 
surface, except the near vicinity of the poles, remains stationary in the laboratory frame. The 
swimming movement is similar to the movement of tank treads, which propel the tank without 
the links dragging on the ground. Therefore, the needle-shaped treadmill (e = d/l <C 1) is 
self-propelled throughout almost quiescent fluid yielding very low viscous dissipation and high 
swimming efficiency. One can not make treadmilling completely non-dissipative because there 
is always some remanent dissipation associated with the motion of the front and rear ends. It 
was demonstrated that the optimal ratio of power invested in dragging and swimming scales 
like (elog l/e)~ 2 and can be made arbitrarily large. 

In this works we consider a closely related mode of locomotion inspired by PurcelFs idea 
PQ and powered by surface tank-treading: "invaginating" torus (FigQ]A.) and two counter- 
rotating disks (similar to FigQjB). This technique, like the treadmilling, involves no change of 
shape, but, opposite to treadmilling, does preserve surface material. In this paper, we will 
address the hydrodynamic model of propulsion, swimming velocity, hydrodynamic efficiency, 
and other issues for three- and two-dimensional swimmer, respectively, We will also propose 
and numerically test a potential design of the microswimmer adopting this propulsion technique 
and composed of rotating microspheres. 



2 Hydrodynamic model 

General solution- In order to calculate the motion of the doughnut-shaped swimmer we intro- 
duce toroidal coordinates [27J (£, rf) in the meridian plane (z, r) via the transformations 

sin 77 sinh£ 

z = c — , r = c — , (1 

cosn £ — cos 77 cosn £ — cos 77 

where 0<£<oo,0<77<27r and c > 0. It follows from (pQ) that the curves £ = const form a 

family of nonintersecting coaxial circles with the centers in the plane z = 0; the typical circle 

£ = £0 has its center at a distance b = c coth£o from the origin and has a radius of a = ccsch£o 

(see Figj2]). Upon rotation about z-axis these circles generate an eccentric family of toruses. The 

toroidal coordinates (£,??,</?) form a right-handed system of orthogonal , curvilinear coordinates 

with metrical coefficients h^ = h, q = h= cosh ^ cosr ' , \ = cos c h s ^ T? . 

We assume that a steady axially symmetric creeping flow (Re <C 1) has been established 

around the swimmer as a result of the tangential surface tank-treading (or "smoking ring") 

motion with a uniform far-field velocity U in the positive z-direction (i.e. in the lab frame the 

swimmer is propelled with velocity U in the negative z-direction) . The flow is governed by the 

Stokes and continuity equations, 

Av = /jgradp, divt> = 0, (2) 
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Figure 2: Toroidal coordinates in a meridian plane; dashed line corresponds to a cross-section 
of a swimmer with b/a = cosh£o = 2, while the torus is obtained upon rotation around z-axis. 

respectively, with the boundary condition at the surface 

v = 14(77) e v at £ = £ • (3) 

For the Stokes flow with axial symmetry the flow problem outlined in (J2H3|) is reduced to 
determination of a stream function ip which satisfies in the region of the flow the equation 
L 2 _ x ijj = E 2 (E 2 ip) = 0, where 

a 2 d 2 kd_ 
oz z or z r or 
The velocity components are readily obtained from ip: 

hdip hdifj 
^ = "r^' ^ = r0£ (4) 

Following [28], where the closely related problem of the uniform flow past an "immobile" 
torus was considered, we seek a solution which is periodic in r\ (of period 2n) and satisfies the 
boundary conditions ([3]) at £ = £0 It follows from and @ that 

dip r 

ip = X, -q7 = u Vl) at £ = & • (5) 
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and ip — > — ^[/r 2 as p — > oo, where p = \/r 2 + z 1 is the distance from the origin to the point in 
space. The constant x should be found as a part of the solution. We follow [28] and write 



-Ur 2 + tp' + xi>i 



(6) 



where ip 1 and ipi give rise to zero velocity at infinity and satisfy ([SJ on £ = £o via the conditions 

(7) 



,. 1 TT o dip' TT dr r 
* =~ 2 Ur , _ = l/r^ + -«fo) 



^1 = 1, 



0. 



(8) 



Note that for "immobile" surface, i.e. u(r/) = 0, the problem is identical to that considered 
in [28] (glider problem). The solution for ip 1 can be constructed from ip 1 and ip 5 = r~ 2 ip , 
where ip k represent a generalized axially symmetric potential corresponding to the solution of 
L k ip k = [23 M, and yields 



^' = r 2 (s - t) 1 ' 2 A n sP^ 1 _{s) + B n P n ^i{s) cosnr ] = r 2 (s-t) 1 / 2 J2W n {s)cosnr ] , (9) 

n=0 ^ J ^ n=0 

where s = cosh^, t = cost/, ^ indicates that the term for n = is to be multiplied by 1/2 and 
^_l/ 2 ( s ) = ^-fn-i/2( s )) with P n being the Legendre function of the first kind 
Substituting ([9]) into the first equation in ([7]) and using the identity 



2\/2 00 ' 

(so ~ t)~ 1/2 = Q n -i ( s o) cosnrj , 



7T 



n=0 



leads to the condition analogous to Eq. (3.15) in [28] 



V2 



A n s P' i (s ) + B n P n _ i (s ) = Q n _i (s ). 



7T 



(10) 



(11) 



Here Q n stands for the Legendre function of the second kind [30 J and so = cosh£o- 

Substitution of © into the second equation in (J7J) and using (fTTj) yields after some algebra 



^ dW n (s Q ) y/2 y—i. , 2 ufa) 

> cos nr/ = (7 > Q i (so) cosnr/ + cscn £o 

z — ' dsn ^ n " 



n=0 



7T 



n=0 



(so-*) 1 / 2 ' 



(12) 



In the particular case of constant tank-treading speed, i.e. u(rf) = au; = v s = Consfl the term 
(s - ty 1 / 2 on the RHS of (02j) can be expanded via ([101) to y ield 

d 



A;- 



ds 



soP'_i(so) 



+ B n P' J (s ) = *7 — i (s ) + — « s csch 2 to Q n _i (so) • (13) 

n 2 TT n 2 7T 2 



For we use the same representation as for ip' in ([9]) with A n and P n replaced with C n /c 2 
and D n /c 2 , respectively, 

9 OO / r 



^1 



n=0 



cos nrj . 



(14) 



the surface rotates clockwise (with respect to (^9-axis) with the angular velocity uj — — uj 



5 



The boundary conditions ([8]) for tpi are identical to those corresponding to the flow past a 
immobile torus, and, therefore, the analysis is identical to that in [28], pp. 85-86. Substituting 
(Q3D into the first equation in ([S]) leads after some algebra to, 

Cn so P' n _ i {so) + D n P n , (s ) = Q" 2 _ ! ( So ) , (15) 

2 2 ixy Z 2 

where Q~ 2 i(so) is the Legendre function of the 2nd kind and degree —2 [30j . The second 

"~2 

condition for ip\ in ([8]) is treated in much the same as way as the analogous condition for ip' 
leading after some algebra to 



C n j — 
dsn 



so PL i(a ) 



3 d 



+ A* i (so) = —s-r-Q^t i (*>) • (16) 
n 2 7rV2ds n ~2 

Therefore, (jlll) . (|13() . (115D and (|16p can be solved for the unknown coefficients A n , £> n , C n and 
D n for all n = 0, 1,2, . . . in terms of the unknown propulsion velocity U and the prescribed 
surface tank-treading speed v s . 

The value of the stream function at the boundary of the swimmer, is found from single- 
valuedeness of pressure p upon integrating it over a closed contour (see [28] for details) 

dp / 1 d 



/^dr=/*f (£V) dr = 
Jc or Jc r on 



where C is any closed contour in the meridian plane, which encloses the body profile, and n 
and r are the unit normal and tangent vector to the curve, respectiveljQ. Integration over the 
contour composed of the segment of the z-axis —R<z<R and the semi-circle p = R that 
joins its end points (±J?, 0) and letting R — > oo yields after some algebra 



X = -c 2 



^An/^Cn). (17) 



vn=0 / n=0 



The swimming speed U should be determined from the requirement of force-free swimming. 
The net force exerted by the viscous liquid in Stokes flows (there is only one component F z due 
to axial symmetry) can be found from the asymptotics of the far-field as [27] 



where V'oo = —l/2r 2 U denotes the stream function corresponding to the uniform flow at infinity 

/ \ 1/2 

(in the frame fixed with the swimmer) and p = v r 2 + z 2 = c I ^£ J . Since p = oo correspond 
to the "point" {£ = 0, r\ = 0}, the above expression for (dimensional) force can be evaluated as 



c 



V2 



OO / 

y 

8irp 2 ^ 

^ n=0 



A n + f 2 C n ) (n 2 -\)+2( Bn+ f 2 D r 



(18) 



where we have made use of identities P'_i(l) = l/2(n 2 - 1/4), = 1 



The sense of t is such that the system of (n, t, e v ) is right-handed in this order 
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Due to the linearity of (|llll3p the solution ip' can be decomposed as ip' = U ip', g \ + v s ip( p y 
The subscript 'g' stands for "glider" : a swimmer with an immobile ( "frozen" ) surface (v s = 0) 
being dragged through the viscous liquid with the velocity {7 = 1; and the subscript 'p' stands 
for "pump": the surface of the torus undergoing continuous tank-treading motion with the 
prescribed tangential velocity v s = 1, while a finite force applied to the torus holds it in place, 
i.e. {7 = 0. Therefore, 

(A n , B n ) = U (A%),BM) + v s (A%\bM) , (19) 
and (|17p can be re- written as 



X = -c 2 



U £ A%) + v s Y, 4 P) = -c 2 («{/ + 0v s ) , (20) 

\ n=0 n=0 / / n=0 

where and can be found by solving llllj) and (|13|) for the "glider" problem (with 
{7 = 1, v s = 0) and the "pump" problem (U = 0, v s = 1), respectively, and C n , D n from 
the complementary problem (fT5Hl6|) of ipi. Further, substituting (fT9|) and (pOj) into (fT8|) and 
equating the net force exerted on the freely suspended swimmer to zero yields the swimming 
speed U as a linear function of the boundary velocity v s : 



U 



y 


'(A$-p Cn ) 


( n2 ~ V 


+ 2 


^i p) - (3D n ) 




T' 


\A^-aC n ) 


i n2 - i) 


+ 2 


[B { n 9) -aD n ] 





x v s . (21) 



With the swimming speed at hand, we can compute the value of the stream function at the 
swimmer surface, x> from (1200 which completes the solution of the self-propulsion problem. 

The expressions for the propulsion speed (|2ip and the viscous force (|18j) can be simplified 
by noting that 

(n 2 - C„, + 22? n = | , (22) 

that can be verified directly using the closed- form expressions for the coefficients^! corresponding 
to the glider problem derived in [28J. 



3 Toroidal glider 

We numerically solve the infinite linear system of equations (truncated at some n max = L) 
corresponding to both, the "glider" and the "pump" problems. Although explicit solution for 
the coefficient A n ~D n can be derived (as in |28j), the resulting expressions are cumbersome, 
and involve integrals like <2"_!/ 2 ( s ) -fn-i/2( s )ds that require numerical integration, we have 
chosen therefore, to use here the numerical solution. The truncation level L corresponds to the 
accuracy of 10~ 6 and varies between L = 45 for b/a = 1.02 and L < 8 for b/a > 2. 

3 Note that ty is missing in denominator of expressions for the coefficients C n and D n appearing in (3.32) and 
(3.33) of Pell and Payne's paper [28] 
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Figure 3: The dimensionless force exerted on the "immobile" torus by the flow uniform at 
infinity vs. the aspect ratio sq = b/a: current theory (solid line); slender-body theory |34j 
(dashed line); numerical results from Table 1 in [32J (•). 



To verify the numerical scheme we first addressed the "glider" case, which is equivalent 
to the well-known problem of the axisymmetric flow past a rigid (i.e. "immobile") torus that 
was studied rigorously [281 Ell Ell E3] and approximately (in the spirit of slender-body theory) 
[34\ I35j. The dimensionless drag force exerted on the torus (scaled with 6irfiaU) calculated as 
a function of aspect ratio sq = b/a is in excellent agreement with earlier results reported in [32] 
for the range 1.3 < sq < 4 and with [31] for the case of closed torus (so = 1). For large values 
of so, the asymptotic form of the force can be obtained following [32] by taking asymptotic 
forms for P n _i/2(so) and Qn-i/2( s o)i solving (fTl]) - (fl"3|) for An and Bn and substituting them 
into (fT8|) together with the identity ([22]) . Alternatively, using slender-body approximation [34J, 
valid at b/a 3> 1, the dimensionless force exerted on a thin torus reads 

F 4 ™° (23) 



6irfiaU 3(log8s + 1/2) 

It can be seen from Figj3] that the expression in (|23p is very accurate even at moderate elon- 
gation, and deviation from the exact result is less than 2% already at sq = 4. As expected, 
in Stokes flows the drag force is controlled by the largest particle dimension and grows (sub- 
linearly) with so as F ~ so/logso similarly to the needle-like slender body. The typical stream- 
line pattern is depicted in Fig 2] for sq = 2. 



4 Propulsion velocity of toroidal swimmer 

The dimensionless propulsion velocity U/v s of the toroidal swimmer powered by the constant 
angular velocity of the surface is calculated via (f2T]) as described in Secj2]and plotted vs. b/a 
in FigE] (solid line). The swimmer is propelled against the direction of its outer surface motion 
(i.e. U > 0) and the scaled propulsion velocity riches maximum, U /v s ~ 0.665 for a closed torus 
(b/a — > 1) and decays as slenderness (b/a) increases. The asymptotic form of the propulsion 
velocity at b/a S> 1 can be derived by considering flow-field produced by line distribution of 
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Figure 4: The streamline pattern corresponding to the flow around toroidal glider in the merid- 
ian plane (the upper-half). The dashed line is the dividing streamline with ip = \- 



rotlets v ro t(x;s) = T 



dC(s) 



x (x — c{s))/(x — c(s)) 3 placed along the torus centerline c(s) with 



ds 



arclength parameter s and density V = l/2u>a 2 [35J. Integrating v ro t over the length of the 
circle in the limit of large sq ° ne obtains the net propulsion velocity [36] 




The agreement between the current theory and the asymptotic result (I24p (dashed line in 
Figj5|) is very good and it can readily be seen that the slender-body approximation provides 
an accurate estimate of the propulsion speed (within 1% accuracy) starting from so ~ 6. The 
scaled propulsion speed decays slowly ~ logso/so for thin ring-like swimmers; for example, 
the swimmer with b/a ~ 20 (e.g. twirling circular dsDNA ring [36]) would translate with the 
speed of ~11.4% of the tank-treading speed. The typical streamline pattern corresponding to 
swimming powered by constant speed tank-treading at b/a = 2 is depicted in Fig[6l 

Stokes flows are rather counter-intuitive: the common intuition suggests that the swimmer 
should be propelled in the opposite direction as the hydrodynamic resistance to the flow inside 
the hole region is higher than outside. This expectation is invalid: the surface in the hole is 
undergoing tank-trading in the same direction with the same velocity everywhere dragging the 
fluid with it through the hole with relatively little resistance, while the outer portion of the 
surface is experiencing higher resistance as it is pushing against the stationary liquid at infinity 

Two interesting questions are (i) whether the toroidal swimmer can perform any better in 
terms of the propulsion speed and (ii) what would be the propulsion speed for a biomimetic 
toroidal swimmer covered by highly deformable but almost inextensible membrane (biological 





(24) 




Figure 5: The scaled propulsion velocity of the toroidal swimmer, U/v s , powered by constant 
velocity tank-treading 71(77) = v s = auJ v = Const vs. the aspect ratio sq = b/a: current theory 
(solid line); asymptotic result of [36] (dashed line); swimmer made of closely packed (along the 
centerline) spheres rotating with constant angular velocity u (•). 

membranes composed of lipid bilayer have a large modulus of dilation, i.e. they behave like 
two-dimensional nearly incompressible media). Since the swimmer moves in the direction of 
the inner surface motion, i.e. it is propelled forward by the outer ("working") portion of the 
surface, one may expect that reduction/increase of the tank-treading velocity in the hole would 
not alter the locomotion speed-wise. In the case of a biomimetic surfac, the inextensibility 
constrain imply that the smaller inner portion of the surface should move faster than the 
larger outer surface. In this case, the surface velocity u should satisfy V s • u = 0, where 
V s = (J — nn) ■ V is the surface gradient operator. This condition in the particular case of an 
axisymmetric membrane becomes |37| 

_ du 1 , , 

V s -u = t ■ — + - u • e r = , (25 

os r 

where r = e v is a unit tangent vector to the surface in a meridian plane, and e r is the unit 
vector in the direction of r. Obviously, the former case of constant tank-treading speed ^ = 
(u = v s e v , v s = Const) is inapplicable for inextensible surfaces. Substituting u = u(rj)e v and 
using the vector identity e v = h(^e r + f^e^) in (|25j) we arrive at the 1st order differential 
equation for u{rj) 

_ du 1 du sin 77 

V s • u = h - — | — u(e r ■ e„) = h u = . (26) 

drj r drj c 

The solution of this equation is straightforward and it defines 77(77) is a unique manner up to a 
multiplicative constant velocity v* that controls the magnitude of the tank-treading: 

uiv) = v s ( cosh Co ~ cos 77) = v s S ° - . (27) 

so + 1 

We chose the constant v* equal to v s (cosh£o + I) -1 so the scaled tank-treading velocity u/v s 
riches a maximum value of 1 at r\ = for any sq, i.e. the surface velocity is scaled with the 
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Figure 6: Flow around the toroidal swimmer driven by constant-velocity-tank-treading motion 
of the surface. The dashed lines denote the closed streamlines in the near vicinity of the 
swimmer's surface. 

maximal boundary velocity v s which for the incompressible membrane is attained at the inner 
surface of the torus. 

The velocity distribution (|27p can also be obtained by much simpler arguments. The con- 
servation of the incompressible surface requires that the 2d flux § c rd(p = Const, where the C 
is any closed circle at the torus surface at £ = £n corresponding to some fixed value of ry. Since 
the surface velocity v = u(rj) is constant along C, the integration leads to 2ttu {rj)r = Const, 
and thus u(rf) ~ 1/r and since at the swimmer surface r = csinh£n/(so — t) we readily arrive 



The typical tank-treading velocity distribution (|27|) corresponding to incompressible surface 
is presented in FigJ7]a for so = 1.1, 2 and 5 vs. the polar angle 9. 9 is measured from r-axis 
about the center of the swimmer cross-section in the upper half-meridian plane, (r = b, z = 0) 
(see Figj2]), so that cos6> = (tso + l)/(so + 0- ^ ^ s readily seen that for a nearly closed torus 
with so = 1-1 the tank-treading velocity of the outer surface at r/ ~ 0, 2tt almost drops to zero. 
Substitution of u(rj) from (|27p into (|12p and some algebraic manipulations yield the infinite set 
of equations analog ous to (|13p for the coefficients A n and B n 



where a n = or |, and (3 n = 1 or ^ according as n = or a positive integer, respectively. 



at d27J). 




+Vs — r -Ct n Q n _Z (S ) + S Q n _I (s ) - PnQ n+ l (So) 

7T Sn + 1 L 2 71 2 Tl+ 2 



(28) 
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u/v s u/v s 




Figure 7: Scaled tank-treading velocity u/v s vs. the the polar angle 9 for sq = 1.1 (red); so = 2 
(yellow); sq = 5 (blue), (a) incompressible surface with velocity distribution (|27|): (6) extensible 
elastic surface with velocity distribution (|291 



We solve the modified pump problem (jlip . f|28[) in the same way as described in Secj2] (the 
glider problem remains unaltered regardless of the particular tank-treading velocity) and plot 
the resulting scaled swimming speed U/v s vs. sq in FigJSJ (yellow curve). The comparison with 
the propulsion speed of the analogous swimmer driven by constant-tank-treading motion (the 
red curve in FigJHJ) suggests that the microrobot propelled by twirling inextensible surface is 
a rather lousy swimmer. Note that re-scaling the surface velocity in (|27p with the speed of 
the outer surface, i.e. writing v* = v s (cosh£o — 1) , yields propulsion speeds U/v s close to 
those corresponding to constant-twirling-velocity swimming (compare the dashed yellow curve 
and the red curve in FigJH]). This supports the hypothesis that the motion of the outer surface 
is in control of propulsion, while tank-treading in the hole region affects only the swimming 
efficiency (see the next section). 

Now, one may also think of a swimmer propelled by tank-treading of highly extensible 
(metamorphic, [2]) elastic membrane, such that the outer surface moves faster than the inner 
surface. Let us assume a simple form of the velocity distribution satisfying this condition 

vt Sq — 1 , . 

u fa = ~^T^ = v s , 29 

cosn £o — cos 77 so — t 

where the non-dimensionalization is with respect to the maximum speed v s occuring in this 
case at 77 = on the outer surface. The typical velocity distribution (|2"9"j) is presented in FigJ7]6 
for various elongations: so = 1.1, 2 and 5. Substituting u(rj) from (|29p into p2|) leads after 
some algebra to 



A ± 

ds 



so PL i(s ) 



+ B n 1 (so) = U — Q' n _i (so) -v s — csch 2 £o Q' n _ 1 (s ) . (30) 

n 2 7T " 2 7T "2 



The resulting propulsion speed of the torus powered by twirling of the extensible surface (129]) is 
provided in FigJH] (blue curve). The scaled propulsion speed is similar to the other two results, 
represented by the red and the dashed yellow curves, respectively, in FiglSJ) for nearly all values 
of b/a. 
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Figure 8: Dimensionless swimming speeds U/v s (scaled with the maximum tank-treading veloc- 
ity of the membrane v s ) vs. the aspect ratio sq = b/a: constant tank-treading speed u(rj) = v s 
(red); incompressible surface with velocity distribution (|27p (yellow); extensible surface with 
tank-treading velocity distribution (j29j) (blue). The dashed yellow curve is corresponds to 
twirling of inextensible membrane re-scaled with the velocity of the outer boundary. 

5 Swimming efficiency 

The results in the previous section showed that the doughnut-shaped swimmer can propel itself 
with a speed of 66% of the twirling velocity of its surface. Now let's determine how efficient 
such swimmer iqj 

The power expended done by an arbitrary shaped organism and dissipated by viscosity in 
the viscous fluid is 

V = - I (a-n)-udS, (31) 
Js 

Eq.flSH) can be written for the axisymmetric tank-treading of torus as — Is 

dS, where 

in toroidal coordinates is given by 

a ^ = h {Cc\ + ^) +^K sinh £ + ^sin7 ? ) • 

Using the fact that for swimmer propelled by purely tangential motion of its surface = 
du^/dn = 0, and applying (JIJ) to determine du n /d^ from tp we integrate a^u^ in (f3T1) numeri- 
cally over 7] from to 2tt with dS = 2irc 2 sinh^oAso — t) 2 ^V- 
The Lighthill's swimming efficiency 5 is defined via 

5 = —p^- , (32) 

where V is the rate of viscous dissipation in swimming with velocity U, and the expression 
in the denominator is the power expended by dragging the glider at velocity U upon action 
of an external force F^ 9 \ For an axisymmetric swimmer F^ 9 '-U = IZfuU 2 , where IZfu 



4 By efficiency we mean only the hydrodynamic efficiency; this does not include the work related to other 
forces, e.g. elastic forces within membrane 
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Figure 9: Swimming efficiency 5 vs. the aspect ratio so = b/a: constant angular velocity u{rf) = 
v s = cllo (red); incompressible surface with velocity distribution (j27[) (yellow); metamorphic 
surface with velocity distribution (|29p (blue); the dashed line correspond to the asymptotic 
result for a "slender torus" ([35]); a swimmer made of closely-packed rotating balls (•) as in 
FigflUJ 

is the appropriate hydrodynamic resistance, this definition reduces to 5 = IZfu U 2 /V. The 
efficiency 5 is dimensionless and compares the particular swimming technique with dragging. 
The higher 5 the more efficient the swimmer is. Using results from previous sections addressing 
the force exerted on the glider and the swimmer's propulsion velocity, we calculate 5 for each 
one of the three modes: constant surface rotation velocity u{rj) = v s = ua, tank-treading of 
incompressible membrane (|27p and twirling of extensible surface (j29|) . The resulting dependence 
5(sq) is provided in FigJ3 As expected, the propulsion powered by the latter mode, i.e. tank- 
treading by extensible surface with nonuniform velocity is superior over the other two cases. In 
the former case, 5 ~ 24% is achieved for a nearly closed torus so — ► 1. The constant twirling 
velocity swimming (the red curve in Figj9|) yields maximal 5 ~ 13% for so = 1.15, which is still 
better than other known propulsion techniques. The swimming efficiency of the torus powered 
by tank-treading of inextensible membrane (yellow curve in Figj9|) is lower, as expected due to 
extra power invested in high-speed tank-treading of the inner surface in the hole, and has a 
peak 5 ~ 3.2% at sq = 1.6. 

Swimming by extensible surface tank-treading is superior over the rotating helical flagellum 
[9], beating flexible filament [6], the Percell's "three- link swimmer" |20tl21|. locomotion by virtue 
of shape strokes [23, 24J and others. For comparison, swimming efficiency of undulating rod 
propagating plane waves of lateral displacement, has an optimum efficiency of 5 = 7.36%, while 
the most efficient swimmer made of two linked counter-rotating helixes of opposite handedness 
has a peak efficiency of 5 = 8.58% [5]. Note that these values correspond to infinite unloaded 
configurations. The swimmer propelled by undulating finite flexible filament attached to a 
passive cargo, as in experiments in [15] is expected to do much worse than 5 = 7.36%; the 
experimentally determined hydrodynamic efficiency of bacterium e. coli powered by rotating 
helical flagellum is ~ 2% |38| . In our case, a useful cargo can be enclosed within the tank- 
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treader without reducing the efficiency. The efficiency of spherical squirmers self-propelled 
by propagating surface waves along their surface (the mathematical model of cianobacteria 
) has the upper bound 5 < |, while numerically calculated values of 5 are usually less 



than 2% [14] : the maximum efficiency reached by the optimized three-link swimmer is only 
1.3% [21j. Although tank-treading cannot, of course, beat the nearly frictionless treadmilling 
technique |26j . where 5 was shown to grow unbounded for a slender needle-like swimmer (as the 
backward motion of the inner surface and possible energy dissipation associated with it were 
not considered), the former mechanism is more attractive from practical point of view. 

The asymptotic behavior of 5 for slender torus, sq 3> 1, can be readily evaluated. To this 

2 

end consider first the flow-field around a straight rod rotating about its axis, v = ^r-e^. Sub- 
stituting (cr-n)-v\ s = (p$*r — ^r) v v = —2au) 2 into (|3"Tj) and multiplying by the surface area 
of the torus, equal to Air 2 ab we immediately arrive at the estimate V ~ 8tt 2 fibv 2 . This estimate 
is, actually, a lower bound for V . Re-writing the dissipation integral as V = 2fi j v E : E dV, 
where E is the rate-of-strain tensor and expressing the product E:E as J2 CiCi + ^■{^i v j)ipj v i)i 
where £ = curli> denotes vorticity, allows expressing V for microswimmers self-propelled by 
purely tangential motional u as |14| 

V = n [ C 2 dV + 2n [ u 2 k s dS . (33) 
Jv Js 

Here V is the fluid volume surrounding the swimmer and k s = —(dr/ds)- n is the curvature 
of the surface along the direction of the flow [14] that is in case of a torus equal to 1/a. Since 
the creeping flow around infinite straight cylinder rotating around its axis is irrotational (i.e. 
C = 0), the volume integral in (|33|) vanishes as b/a — > oo. At finite b/a we have, however, only 
an inequality 

V>2^i[ u 2 K s dS = 8tt 2 fibv 2 . (34) 
Js 

Comparing with the asymptotic form for the viscous force on a glider (|23p and with the propul- 
sion speed (f24"|) . one can readily estimate 5 from above as 

{< (log 8s „-i) 2 , 
- 4 5 g(log8»„ + i) 14 

This asymptotic result is presented in Figj9] as a dashed line and is in a very close agreement 
with the exact result already for moderate so: the relative error between the asymptotic and 
the exact result is ~ 1% for sq = 5. 



6 Microswimmer made of rotating spheres 

In this section we propose a potential design for the swimmer propelled by the surface rotation. 
Imagine the necklace-like ring of radius b in the xy plane made of iV closely packed (separated 

5 In the laboratory frame the velocity at the surface is a superposition of the translational velocity U and 
purely tangential motions u 
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Figure 10: Illustration of necklace-like swimmers made of rotating spheres. 

by the distance d = 0.05a) spheres of radii a (see Fig JIOp . In a cylindrical polar coordinate 
system (z,r,(p), each sphere rotates at the constant velocity £2 = — fie^ , which in the absence 
of external forces causes the necklace to swim, as described in the previous sections, along the 
normal to the plane of the ring in the negative z direction. This concept may serve as a basis 
for potential experimental design of a micro (nano) robot based on Purcell's idea. 

Calculations of the Stokes flow involving collection of spheres are rather standard today. 
We use the Multiple Expansion (ME) approach [40J and construct the solution of the Stokes 
equations in terms of Lamb's spherical harmonics expansion. The no-slip conditions at the 
surface of all particles are enforced rigorously via the use of direct transformation between solid 
spherical harmonics centered at origins of different spheres. The ME method yields a system of 
3xiVxLx(L+2) linear equations for the expansion coefficients and the accuracy of calculations 
is controlled by the number of spherical harmonics (i.e. truncation level), L, retained in the 
series. In some situations, when the rigid particles in near contact approach each other, the 
hydrodynamic interaction involve strong lubrication forces (~ 1/d singularity in the viscous 
forces for "squeezing" flow between two spheres) and the ME method requires a large number 
of spherical harmonics due to slow convergence. The "near-field" (closed-form) lubrication 
"patches" based on pair-wise hydrodynamic interaction [39] can readily be introduced to greatly 
improve convergence. In our case, the hydrodynamic interaction between neighboring spheres 
is weaker as the distances between the spheres are kept fixed, and only involve shearing of the 
fluid in the thin gap produced by their rotation (the viscous force ~ log (1/d) for "sliding" 
motion between two spheres). L < 7 was sufficient for all configuration with N < 28 spheres 
to achieve an accuracy of less than 1%. 

The resulting propulsion speed scaled with aSl is depicted in FigJJ] (• symbols) for swimmers 
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composed of from N = 2 to 28 rotating spheres. The maximum propulsion velocity of U/aQ — 
0.37 is obtained for "quadrocycle" (i.e. for N = 4). Though the swimmer made of spheres 
moves slower than the rotating torus (see FigE|), the propulsion speed is yet considerable. 

For necklace-like swimmers made of ./V identical microspheres as in FigllOlthe rate-of- viscous- 
dissipation is calculated as V = J2j T^j = NV sv h where V sp h corresponds to dissipation rate due 
to a single sphere in the necklace, jth sphere in the necklace undergoes a rigid-body motion 
according to Va = U + fl X rj with rj is the radius vector with origin at the particle center, U 
is the swimming speed and fl is the angular velocity, and thus 

V sph = - (cr • n) ■ v dS = -U-F sph - fl-L sph , (36) 
JSj 

where L S ph = 

Is r 3 x ( (T - n ) dS is a hydrodynamic torque exerted on each (jth) sphere of 
the swimmer. Since each sphere in the necklace is force-free, F sp h = 0, and the swimmer is 
torque-free, J2j(Lj + Rj x Fj) = 0, (Rj is the vector connecting the swimmer origin and the 
center of the jth sphere) the propulsion efficiency (|32p can be determined from 



F { s 9 pl U 



■'sph 



Kfl}U 2 _ a 2 Tlfl} IV 

Tjsph 2 TD s P h 



(37) 



where TZ S p^ and are the appropriate hydrodynamic resistances functions (configuration 

dependent) for a single sphere in the necklace. These functions can be readily calculated 
numerically via the ME method: TZ S p^ ~ \xa is equal to the viscous force exerted on a single 
sphere in the glider dragged through viscous liquid with velocity {7 = 1, while ~ /^a 3 is 
equal to the torque L exerted on a single sphere in the freely suspended swimmer with all sphere 
rotating with the angular velocity — £le v . The results are depicted in Figj9]for N = 2-12, while 
the maximum, 5 ~ 4.95% is obtained, again, for "quadrocycle". 



7 A two dimensional swimmer 

A natural two-dimensional analog of the toroidal swimmer consists of two discs counter-rotating 
at the same angular velocity u. As we show below this has a simple closed form solution. A 
realistic swimmer based on this description may be constructed by attaching two counter- 
rotating cylinders. 

The two-dimensional Stokes equations are conveniently dealt with by employing the complex 
notations v = v x + iv y , z = x + iy and d = ^(d x — id y ). In these notations the equations take 
the form 2fiddv = dp, Ke(dv) = and the general solution is given by v = f + g — zf, p = 
— 4/iRe(/') with /, g being arbitrary holomorphic functions (e.g. [41J). 

Finding /, g explicitly is, in general, a difficult task. In order to solve the problem at hand of 
two counter-rotating discs of radius a with their centers separated by a distance 2b we employed 
conformal mapping techniques. The details of these lengthy calculations are presented in the 
Appendix together with an explanation on how to extend it to a more general (nonuniform) 
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boundary velocity on the discs. The final result, however, turns out to be fairly simple and is 
given by: 

a?v ( z 2 — 6 2 + a 2 \ 

f{Z) = 46(6 2 -fl 2 ) ^2 + 6 2_ fl 2 ) > ( 3§ ) 

[ a(26 2 -a 2 ) q(46 2 - a 2 ) a^ 2 - a 2 ) 1 

Us \ 46(6 2 - a 2 ) + 26(z 2 + 6 2 - a 2 ) b(z 2 + 6 2 - a 2 ) 2 J ' 1 j 

where u s = au is the boundary velocity. Indeed, substituting z = ±i& + ae 1 ^ one readily finds 
v = iri^e 1 ^. Since /, g are regular everywhere outside the discs (including at infinity) one 
concludes that this is the correct flow solution. 

At infinity we have v(oo) = —-^v s i.e. the fluid flows in the direction where the outer 
portion of the discs moves as in the 3D case (here the negative ^-direction). The swimming 
speed, in the fluid rest frame, will then be 

U = +^v s e x . (40) 

The element of force acting on an infinitesimal segment may always be expressed as dF = 
\pdz + (2i/i <9v) dz = 2i/xd(u — 2f). Since / in our solution is single- valued (no log term) 
one immediately concludes that the two discs exert no force on each other: § dF = 0. This 
conclusion can also be derived from symmetry arguments by considering the effect of reversing 
the rotation direction of the discs. 

It is clear that to maintain the rotation, the discs must exert on each other a non-zero 
torque. Straightforward calculation yields the torque on the disk 

L dlsk = ImfzdF = ^av s . 

The dissipation occurring on each disc is given by 

'Pdisc = —Re j> V dF = ATTflVg = Llo . 

The total power on the two discs is V = IVdisc- We recall that in 2D the dragging problem admit 
no regular solution within the Stokes approximation^. Thus defining the swimming efficiency 
as 5 = (F^ 9 ' ■ U)/V makes no sense in the 2D context in which is not defined. This may 
be considered as a mere issue of normalization. A natural measure for the efficiency of the 2D 
swimmer is given by the dimensionless ratio |24} [26] 

_ 4vr/x[/ 2 1 /a x 2 



V 8 V6 



6 This is known as the Stokes paradox and can be resolved by noting that far from the object the quadratic 
term (v ■ V)« ~ (U ■ V)v ~ fiAt! cannot be neglected 
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8 Discussion and conclusions 



In this work we considered a propulsion technique that involves no change in the swimmer's 
shape and is powered by surface "tank-treading". The underlying mechanism of propulsion is 
based on the difference in effective viscous friction to the rotation of "inner" (i.e. in the hole) and 
"outer" torus surface. This idea for potential locomotion at low Reynolds number was originally 
proposed by Purcell pQ more than 30 years ago, but the adequate mathematical analysis (in 
terms of propulsion speed and efficiency) has not been provided until now. Exact solution to the 
Stokes problem for the doughnut-shaped swimmer was obtained here via expansion in toroidal 
(ring) harmonics; For the analogous 2D swimmer the explicit closed form solution was found 
using complex variable techniques. 

It is demonstrated that the doughnut-shaped swimmer outperforms its bio-inspired competi- 
tors - swimmers powered by undulating or rotating flagella, both speed-wise and efficiency- wise. 
The swimmer moves in the direction of the "inner" surface motion: the motion of the outer 
surface yields viscous thrust as it pushes against the motionless liquid far from the swimmer, - 
the opposite from what common intuition would expect in inertia-dominated environment. (This 
intuition might be wrong even when inertial forces takes over: a vortex ring in the inviscid fluid 
is propelled in the same direction as the toroidal swimmer in viscous liquid [42J ) . The swimming 
speed is dominated by the surface velocity on the outer side of the torus, the inner side having 
only small affect. For nearly closed torus propelled by tank-treading of extensible (biomimetic) 
membrane it can get up to 66% of the surface velocity (for constant surface twirling velocity). 
The hydrodynamic efficiency is more sensitive to the swimming mode. It can get up to ~24% 
for nonuniform surface twirling velocity and probably even higher for appropriate modes. As 
the hole radius, b, increases (for a fixed cross-section radius a), the propulsion velocity slowly 
diminishes as ^- ~ lo y , where so = b/a and v s denotes the tank-treading velocity of the 
outer "working" portion of surface. The swimming efficiency of the doughnut-shaped swimmer 
decays as 5 ~ log .f° for sq 3> 1. When we assume that the torus is powered by tank-treading 
of incompressible surface (bilipid membrane), the surface velocity of the inner (smaller) sur- 
face is significantly higher than that of the larger outer surface. Thus, the propulsion speed of 
such swimmer is considerably lower than that of the swimmers with extensible surface and it 
amounts to 19% of the maximum tank-treading speed v s at b/a = 3.2. The optimal hydrody- 
namic efficiency is 5 ~ 3.19% at b/a = 1.6, but it is yet comparable with efficiency of flagellated 
microorganisms |38j . 

Since there is, probably, no simple route of building the doughnut-shaped tank-treader, we 
propose a practical design of the swimmer that mimic tank-treading propulsion mechanism. It 
is composed of N rigid nearly touching spheres of radius a assembled in a circle and rotating 
about the line of centers with common angular speed f2. The resulting hydrodynamic problem 
is resolved numerically using the Multiple Expansion Method. The "quadrocycle" swimmer 
composed of four spheres is optimal in terms of both, speed, U ~ 0.37 a£l, and hydrodynamic 
efficiency, 5 ~ 4.95%. 
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For the two-dimensional analogue of the invaginating torus, i.e. two close counter-rotating 
disks, we obtain an explicit solution for the hydrodynamic problem, the resulting swimming 
speed and hydrodynamic efficiency, being U = and, 5 = ^p. Therefore, as expected, a 
torus is a better swimmer than its two-dimensional equivalent: the propulsion velocity for the 
2D swimmer (|40p and the asymptotic expression for slender torus at sq ^> 1 ()24p differ by a 
multiplicative factor of logso- 
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Appendix 

Here we consider the 2d flow around a pair of discs of radius a centered at (x, y) = (0, ±6) and 
counter-rotating with angular velocity to. Our aim is to resolve the viscous (Stokes) flow around 
the disks. However most of the derivation can be extended to the case where an arbitrary 
velocity w is prescribed as boundary condition on the two circles. We therefore present the 
derivation in a way that makes this extension obvious. We do, however, restrict ourselves to 
flow being symmetric with respect to the x-axis. 

Using complex notations z = x + \y,v = v x + \v y the general solution to the Stokes equations 
is expressible as v = f + g — zf ,p = — 4//Re(/') where /, g are arbitrary analytic functions (not 
necessarily single- valued). Since we assume the flow is symmetric with respect to the x-axis we 
have v(z*) = v(z)* which imply also f(z*) = f(z)*,g(z*) = g(z)* . 

It is convenient to define new parameters a,/3 by b = , a = P 1 2 ^ a a (equivalently 

j3 = Vb 2 — a 2 , a = \{b — Vb 2 — a 2 ) < 1). We also define a new complex coordinate ( by 
z = i/?i±£. The area outside the discs is mapped by this change of variables to the region 
a < |C| < a _1 (in particular z = oo is mapped to C = 1)- Standard bi-polar coordinates may 
be defined by £ + ir/ = log£, however we will not make use of these. In terms of £ the general 
solution takes the form 

v = f + g + l^(l-C) 2 f' 

The functions /(C)>5 f (C) may contain a log£ term in case they are multivalued. But apart 
from this, they are analytic functions in the annulus a < |£| < a -1 and may, therefore, be 
Laurent expanded as /(C) = E~=-oo a n ( n + log(. . .), g(() = En=-oc + log(. . .)• The 
symmetry of the problem implies a_ n = a*,6_ n = 6*. It is convenient to denote by F,G the 
sums of positive powers appearing in the expansions of /, g. They are clearly analytic inside the 
disc |C| < a" 1 . We also denote F(x) = ££°=i<z n = F(x*)*,G(x) = En=iKx n = G(x*)*. 
Thus, the functions /, g can be expressed as: 

/(C) = F(C) + F(l/C)+ia logC (Al) 
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9(0 



G(C) + G(l/C)-ia lo g C + &o. 



(A2) 



The constant 60 was arbitrarily appended to g, while the reflection symmetry implies it must 
be real. The constant ao must also be real in order for the physical flow v to be single valued. 
It may be remarked that the case under study (i.e. two rotating discs) posses also a symmetry 
upon reflection about the y-axis which may be used to show that all the coefficients in the 
expansion of f,g are real and in particular ao = 0. However, for the sake of generality we shall 
not use this symmetry here. 

Next the general solution v(Q should be matched with the prescribed boundary condition. 
Due to reflection symmetry it is enough to match boundary condition on one circle. We, 
therefore, consider in the following only |£| = a and use the relation £* = a 2 /£ to express v as 
a function of £ alone. We find 



F(()+F(1/0 + g(jJ +G^)+b + ia loga 2 + 



(A3) 



+ 11±C 
2 1-C 



2 r 



1-- 



c 




a- 



a- 



This should equal the given boundary condition u?(£). For rotating discs we have w(z) = 
i^f (2 — ia) implying w(Q = v s ( "^-c" ) ' wnere v s = auJ - 

We now operate on both sides of equation (jX3]) with: §\c\ =a ^ t^a£ • ^ ne analyticity of 



\C\=a 2th C- 

F, G, F, G inside the disc |£| < a^ 1 , allows using the residue theorem to compute the integrals 
involving F(Q,G(-^),F'(-^). The integrals involving F(±),G(^),aF'(^) may be calculated 
in a similar fashion by first changing the integration variable to £' = l/£. 
There are two distinct cases to consider. In the case of |£| < 1 one finds: 

dc w(o _ F( „^ , (a-£) 2 (l+a£) 

(\= a 2-k\ ( - a£ 



2a 2 (l - a£) 



F' 



+ i 



a 



Ota® 



(A4) 
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-1— + b + ia log a 2 - ^F'(O) + ^ r 

2£ 2 1 — a£ 



-F'(a^) . 



For the rotating discs we have w\ (£ ) 



a-* 



l-a£, Ug - 

In the case of of |£| > 1 the same integral becomes 



-F 



(a-£) 2 (l + a£ ) 5, 
2^ 2 (1 - a£) 



a 



a 



2n2 



a£ 



^( a 2)--F'(0)+i^. 

1 v ; 2 w 2£ 



We find it more convenient here to redefine £ as £ — ► l/£* so that the new £ again ranges 
inside the unit disc. It is also more convenient to consider the complex conjugated expression. 
Then we can write 



w 2 (0 



w(C) 



«/£* 



-F 



£\ ^/ « 1 « + £/-, 



a£) 2 F'(a£) (A5) 
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a — £ 2 2 
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For the counter-rotating discs W2 = 0. More generally one may verify that if w(ae lS ) = 
E n w n e m9 then Wl (g) = En= w n ^ n , w 2 (0 = - En=l w ~nC- Eliminating G from ea (IA4IA5l) 
we arrive at 

W(0 = wtfa) + wa(e/a) = ^(« 2 - F (J,) + glgfSgg) ^^) (A6) 
+(1 - a 2 ) 2 y^F'ia 2 ) + ^^'(a 2 )) + ia (l - a 2 )^ - ReF'(O) + 6 + iao log a 2 

For counter-rotating discs W(£) = "^"^ t> s ■ 

We obtained a non-local ODE for the function F. Note, however, that apart from the 
constant F'(a 2 ) = (^F'{a 2 )^j all occurrences of F contain a multiple of £ as an argument. 
Thus, one may hope that expanding in powers of £ around the origin will give useful relations 
that will allow determining the coefficient of F(x) = a nX n in the Taylor expansion at least 
recursively. The first two relations obtained this way deserve to be presented explicitly: 

W(0) = (1 - a 2 ) 2 F'(a 2 ) + b - ReF'(0) + ia log a 2 - ia (l - a 2 ) 
W(0) = (1-a 2 ) 2 



o?F>{a 2 ) - ^F'(a 2 
cr 



iao«(l — a 2 ) 2 



The first relation can be solved for F'(a 2 ) (thus also giving its complex conjugate F'(a 2 )). 
Substituting this into the second relation and separating the equation to its real and imaginary 
parts allows then to solve for ao, bo . Using this to eliminate F'(a 2 ), F'(a 2 ), ao, bo from Eq. (|A6|) 
then leads to a considerably simpler equation. For the counter-rotating discs the equation takes 
the form 

(a 2 - 0(1 " « 2 £) \f(& 2 ) ~ F{i/a 2 )} + (1 - a 4 )£(l - 0^(0 = «U " <*>) v s f (A7) 



In the more general case the non-homogeneous term on the r.h.s of (|A7|) is replaced by a more 
complicated expression constructed from the prescribed W. 

Substituting F(x) = J2^=i in^™ i n (| AT|) one may recursively solve for all the coefficients 
a n in term of a\. Indeed, this procedure appears to work also when the r.h.s is replaced by 
the expression corresponding to an arbitrary boundary condition w. For the special case of the 
rigid-body rotation with constant angular velocity the procedure leads to particularly simple 
solution for a n 's that allows summation F(x) = Y^ a n,x n - One finds 

t?( \ a3U ° >n " 



(l-a 2 )(l-a 4 ) 1-x 

where C is an arbitrary constant. Since F(Q is (by definition) analytic inside the disc \(\ < a" 1 
the constant must vanish, C = 0. Using the other relations we have it is then a simple matter 
to find 

a(l — a 2 + a 4 ) 



a Q = 0, b 
G{x)~- 



l-a 2 )(l-a 4 ) s 
ax(2 + 2a 4 - 2a 2 + a 2 x) 



2(l-a 2 )(l-a 4 ) 



22 



Using (|A2|) and transforming back to z-coordinate we finally arrive at the Eqs. (|38p (|39j) . 



f( z ) = 



9(z) 




a 3 (b 2 - a 2 ) 
b{z 2 + b 2 - a 2 ) 2 



A comment on the case of arbitrary boundary condition w. the iterative solution of the 
generalized Eq. (|A7p depends on the chosen value of a\. Since the correct value of a\ is a priori 
unknown the obtained F{x) may be different from the correct solution up to a homogeneous 
solution with an arbitrary constant C. This issue may be resolved in two ways. The first 

(and most straightforward) is to notice that C = lim^^oo a n . One may use this to estimate 
C and correct the solution accordingly. The other potential approach relies on the fact that 
although is a solution of the homogeneous equation (|A7p with v s = 0, it is not a solution 

of (|A6j) with W = (i.e in the passage from (IA6D to (IA7D some information was lost.) 
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